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We present a theoretical study of ac charge transport arising from adiabatic temporal variation of zero-field spin splitting 
in a quasi-onedimensional hole system (realized, e.g., in a quantum wire or point contact). As in conduction-electron 
systems, part of the current results from spin-dependent electromotive forces. We find that the magnitude of this current 
contribution is two orders of magnitude larger for holes and exhibits parametric dependences that make it more easily 
accessible experimentally. Our results suggest hole structures to be good candidates for realizing devices where spin 
currents are pumped by time- varying electric fields. 



In conventional electric circuits, flow of charges is induced 
by applied electric -potential differences. As an alternative 
particularly suited for miniature electronic devices, pumping 
of charge by manipulating intrinsic system parameters has at- 
tracted a lot of attention. Several intriguing proposals of this 
kindi^ are based on the tunability 3 of spin splittings induced 
by structural inversion asymmetry (SIA) in semiconductor 
nanostructures. 4 Such exotic effects could be studied exper- 
imentally in materials with intrinsically large spin splitting. 
Generally, p-type structures are favorable 4 ' 5 because valence- 
band states are directly affected by the atomic spin-orbit cou- 
pling, while conduction-band states experience a sizable SIA 
spin splitting only in narrow-gap materials. Also, the larger 
effective mass of holes implies that spin splitting is a larger 
fraction of the Fermi energy at given carrier density than in 
electron systems. 

Here we present a study of charge transport induced by 
time-dependent SIA spin splitting in ;?-type (hole) quantum 
wires or point contacts, taking into account the special prop- 
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FIG. 1. Inset: Setup of a charge-pumping device. Constant side- 
gate voltages are used to adjust the point-contact width w. Time- 
varying top- and/or back-gate voltages modulate the SIA (Rashba) 
spin-splitting strength (measured in terms of a voltage V z ) in the 
quantum well of width d. Main panel: Dispersion relation of the 
lowest quasi- ID subbands in zero magnetic field for aspect ratio 
w/d — 1, using band-structure parameters applicable to GaAs 
[7 = (272 + 37 3 )/(57i) = 0.37] and V z = 0.1 Vr. The volt- 
age scale Vr = 7i7r/i 2 /(2mo |rfi 8 " |) quantifies the magnitude of 
spin-orbit effects in a material's band structure. 



erties of hole states due to their effective spin-3/2 degree of 
freedom. 6 Figure [T] illustrates the basic device setup consid- 
ered in our work. A quasi-one-dimensional (quasi- ID) hole 
system with free propagation direction for holes parallel to the 
x axis is realized by engineering the appropriate quantum con- 
finement in the orthogonal in-plane (y) direction of a p-type 
semiconductor heterostructure that has growth direction par- 
allel to the z axis. SIA spin splitting is modulated in time by 
means of front- and back-gate voltages 3 that generate an elec- 
tric field parallel to the heterostructure-growth ( z) direction. 
As for conduction electrons, 1 the time-dependent SIA spin 
splitting acts as a time-dependent and spin-dependent gauge 
field. The resulting spin-dependent electromotive forces gen- 
erate oppositely directed and spin-polarized partial currents 
from spin-split subbands. To convert the induced spin cur- 
rent into a more easily measurable charge current, an in-plane 
magnetic field is applied perpendicularly to the quasi- ID sys- 
tem (i.e., in y direction). As we will see below, the larger 
magnitude and favorable sample-parameter dependences of 
this current in hole systems make it more easily experimen- 
tally accessible than in n-type structures. 

We use the Luttinger model 7 for the top-most valence band 
in common semiconductor materials. Within the spherical ap- 
proximation, - the corresponding effective-mass Hamiltonian 
for holes reads (we count hole energies as positive) 



H, = 
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Here mo is the electron mass in vacuum, and the parameter 
7 = (272 + 373)/ (571) measures the splitting between heavy- 
hole (HH) and light-hole (LH) bulk-valence-band dispersions 
in terms of the materials-dependent Luttinger parameters.— k 
denotes the hole wave vector, and J is the vector of spin-3/2 
matrices, 4 The Zeeman coupling of hole spin to a magnetic 
field B is described by Hz = 2k/jbB • J, where /ie denotes 
the Bohr magneton and k is the isotropic hole g-factor.— SIA 
gives rise to a coupling between the holes' orbital motion and 
their spin. The corresponding (Rashba) term in the Hamil- 
tonian reads 4 Hr = rf" 8 " (k x £) J. The band-structure 
parameter rf j 8 " = 14.62 eA 2 in GaAs (values for other semi- 
conductor materials can be found in Ref. 4), and S is the ef- 
fective electric field quantifying SIA. 
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(fc 2 ) = n 2 ir 2 /w 2 , with integer m and n. This 



Following the spirit of previous proposals for quantum 
pumps of conduction-band electrons, 1 ' 2 we consider the holes 
to be confined in a quasi- ID system by a quantum- well con- 
finement in z direction and additional in-plane confinement 
V(y). For simplicity, we assume the confining potential in 
both directions to be of hard-wall type, with respective widths 
d and w, and project onto quasi-lD bound states by setting 
k z -> (k z ) = 0, k y -> (k y ) 
and k 2 -> 

yields an effective Hamiltonian H^ mn \k x ) describing HH- 
LH splitting and mixing within the subspace of a single quasi- 
lD bound-state level. As SIA will typically only be sizable 
in the quantum- well growth direction, we have £ — £ z z. 
Defining the abbreviation J± = (J x ± iJ y )/y/2, as well as 
a natural energy scale Eq = ir 2 h 2 ji/(2niod 2 ), magnetic- 
field strength B = j 1 n 2 h/(2Ked 2 ) w 18.9kT/(d[nm]) 2 
for GaAs, the SIA voltage scale V R = ^irh 2 / {2m \rlf v \) 
(sw 5.72 V for GaAs), and V z = £ z d, we find H^ mn \k x ) = 



E [n^ + h$> + n w (k x ) + Wrz(*»)], with 
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Straightforward diagonalization of i/( m ") yields the spin-split 
quasi-lD hole-subband dispersions. See FigureQ] 

HH-LH mixing between different orbital bound-state lev- 
els is neglected within our model. Such a simplified approach 
yields reliable results only for the lowest spin-split subband. 1 1 
Hence, we limit our study to the case where only the two sub- 
bands with m = n = 1 are occupied. Thus four propaga- 
tion channels exist for holes in the device. Their correspond- 
ing Fermi wave numbers, denoted by k[ + 2 for right-movers 

and k[ ^ for left-movers, are straightforwardly found from the 
subband dispersions. 

To model the not necessarily perfect transmission of a hole 
point contact, we introduce a 5-barrier potential Vq&(x) half- 
way between the two leads that make contact to the quasi- ID 
hole system of length L. The scattering matrix S is found by 
adapting standard methods 12 to the case of four-spinors de- 
scribing spin-3/2 hole states. 

Application of (front and/or back-gate) voltages can be used 
to effect a change in the SIA spin splitting 3 and should thus 
make it possible to impose a variation of V z with time r. Un- 
der certain conditions, transport of charge and/or spin will 
be the result of such a time dependence. 1 ' 2 To produce a dc 
charge current in the adiabatic regime, temporal variation of 
at least two parameters is necessary, whereas a time-dependent 
V z by itself only generates ac currents. Studying the latter will 
provide useful insight into the physical mechanism for spin- 
dependent quantum pumping and is also more easily possible 



experimentally. Hence we focus here on ac charge transport 
generated in situations where only SIA spin splitting is time- 
dependent. 

Using techniques introduced in Ref. [13| (see also Ref. 14), 
currents injected in the leads are given in terms of contribu- 
tions due to transmission and reflection coefficients. We label 
different blocks of the scattering matrix by s G {r, r',t, t'} 
and define partial currents i s a*^ 
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Currents entering the left (L) and right (R) leads are then 
given by i L (r) = i t >(r) + i r (r) and i R (r) = i t (r) + v(r). 
We consider the linear combination of left and right currents 
that is related to net charge transport between the two leads: 

Wr) = [i L (r) - *r(t)]/2. 

For our chosen device geometry, ilr is entirely due to spin- 
dependent electromotive forces. These arise because the wave 
function of charge carriers travelling between the leads ac- 
quires a dynamical phase that appears in the transmissions: 

tjj oc expfikjL) and t'^ oc exp(— ifcj L). As the wave 

vectors kjj are functions of V z , they are time-dependent 
when SIA spin splitting is modulated. Time-dependent phases 
of transmission coefficients can be interpreted as electro- 
motive forces that, in the present case, are also subband- 
dependent. This is the origin of the spin-orbit gauge-field con- 
tribution to i t (,) , which scales linearly with L. It turns out that 
all other contributions to zlr cancel, and that the current ilr 
is antisymmetric in B y . The antisymmetry of zlr w.r.t. B y is 
the distinctive feature of charge transport generated by time- 
dependent spin-orbit gauge fields. This is so because spin- 
dependent electromotive forces generate oppositely directed 
partial currents from the two spin-split subbands whose im- 
balance in B y ^ results in a finite transport current ilr. 

Figures [2] and [3] show the dependence of the spin-orbit 
gauge-field generated current «lr on a number of experimen- 
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FIG. 2. Charge current between L(eft) and R(ight) leads generated 
by time-varying spin-orbit gauge fields, plotted as a function of the 
aspect ratio d/w for values of the Fermi energy Ep = 1.0 -Bo (dot- 
dashed curve), 1.1 Eq (dashed curve), 1.2 Eo (solid curve). Other 
parameters are y — 0.37 (the value for GaAs holes), Vo = 0.1 Eod, 
V z = 0.1 V R , B x = 0, and B y = 0.02 B . i = eV z /V R . Currents 
are scaled by the factor d/L. 
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FIG. 3. Dependence of ac charge current on the strength Vo of the S- 
function barrier that models non-unitary transmission of a point con- 
tact (solid curve). The current unit is io = eV z /Vs., and parameters 
are 7 = 0.37 (the value for GaAs holes), V z = 0.1 Vk, E F = 1.2 £0, 
Si = 0, B y = 0.02 Bo, and w/d = 2. For comparison, we also 
show the current (multiplied by a factor of 200 to be visible in the 
plot) generated in a GaAs electron system for the same V z , Ef, and 
B, with band nonparabolicity included (dot-dashed curve) and ne- 
glected (dashed curve). Currents are scaled by the factor d/L. 



tally relevant quantities. We use io = eV z /VR as the unit of 
current and scale current values by d/L to absorb the linear 
dependence on L. 

The magnitude of «lr as a function of the aspect ratio for 
the quasi-lD confinement is shown in Fig. [2] While there is 
a precipitous drop in current as d/w deviates from one, it still 
remains finite in situations where the nature of the propagating 
states becomes more 2D-like. 

Figure [3] elucidates the dependence of zlr on the 5-barrier 
strength and, thus, on the transmission of the quasi- ID system. 
For comparison, we show our result for a hole system along- 
side two curves that are obtained for conduction electrons in 
the same material under the same conditions, with and with- 
out band nonparabolicity included. Two major differences 
between the conduction-electron and hole cases are apparent. 
Firstly, in the same semiconductor material, the maximum of 
zlr for a /?-type structure is two orders of magnitude larger 
than that for an n-type one. Secondly, the dependence on the 
barrier potential is qualitatively different in the two cases. For 
holes, the maximum zlr occurs in the limit Vo — > when the 
quasi- ID system is perfectly transmitting. In a point contact, 
this situation can be realized over a wide range of system pa- 
rameters at the lowest conductance plateau. In contrast, the 
maximum transport current in conduction-electron systems 
appears for a non-unitary transmission value, and reaching it 
in a point contact requires fine-tuning of parameters in the nar- 
row conductance-step region. 

The strong enhancement and different qualitative behavior 
of zlr for holes as compared with electrons can be traced back 
to the strong mixing between HH and LH states in a quasi- ID 
system. The resulting nonparabolicity of hole subbands in- 
creases the sensitivity of Fermi wave vectors and, hence, of 
quantum phases in the transmission amplitude, to time de- 
pendence of SIA spin splitting. Due to their special nature 
and comparable energy scales, HH-LH mixing and HH-LH 



splitting result in qualitatively different nonparabolicity ef- 
fects than those present in conduction-electron systems. To 
illustrate this point, we show in Fig.[3]the transport current that 
is obtained when nonparabolicity in the conduction bandi&ii 
is taken into account 18 (dot-dashed curve) alongside the result 
for the same system with nonparabolicity neglected (dashed 
curve). It is seen that, for conduction electrons, nonparabol- 
icity also enhances the ac charge current generated by spin- 
dependent electromotive forces, but the enhancement is only 
a fraction of the maximum current obtained in the parabolic 
approximation, and the general shape of the Vo -dependence is 
similar in the cases with and without nonparabolicity. In par- 
ticular, the current maximum still occurs at approximately the 
same finite value of barrier strength Vq. This moderate influ- 
ence of nonparabolicity in the conduction-electron case has to 
be contrasted with the drastically different shape of the curve 
obtained for a hole system. 

In conclusion, charge currents flowing in quasi- ID hole 
systems in response to a time-varying SIA spin splitting are 
predicted to have a large magnitude (in GaAs, 200-times 
larger than for conduction electrons) due to strong spin-orbit 
effects in the valence band and because of subband non- 
parabolicity arising from HH-LH mixing. Proportionality to 
system size indicates the purely quantum origin of these cur- 
rents, and their odd symmetry with respect to transverse in- 
plane magnetic fields enables unambiguous extraction of this 
special current component in experiments. 

'M. Governale, F. Taddei, and R. Fazio, Phys. Rev. B 68, 155324 (2003). 
2 A. G . Mal'shuko v, C. S. Tang, C. S. Chu, and K. A. Chao, 

Phys. Rev. B 68, 233307 (2003) Y. Avishai, D. Cohen, and N. Nagaosa, 

Phys. Rev. Le tt. 104, 196 601 (2010)| 

^G. Engels, J. Lange, T. Schapers, and H. Lilth, Phys. Rev. B 55, R1958 
(1997); J. Nitta, T. Akazaki, H. Takayanagi, and T. Enoki, Phys. Rev. Lett. 
78, 1335 (1997); D. Grundler, Phys. Rev. Lett. 84, 6074 (2000). 

4 R. Winkler, Spin-Orbit Coupling Effects in Two-Dimensional Electron and 
Hole Systems (Springer, Berlin, 2003). 

5 R. Danneau, W. R. Clarke, O. Klochan, A. P. Micolich, A. R. Hamil- 
ton, M. Y. Simmons, M. Pepper, and D. A. Ritchie, Appl. Phys. Lett. 
88, 012107 (2006); O. Klochan, W. R. Clarke, R. Danneau, A. P. Mi- 
colich, L. H. Ho, A. R. Hamilton, K. Muraki, and Y. Hirayama, Appl. 
Phys. Lett. 89, 092105 (2006); L. P. Rokhinson, L. N. Pfeiffer, and 
K. W. West, Phys. Rev. Lett. 96, 156602 (20 06); J. Sh abani, J. R. 
Petta, and M. Shayegan, |Appl. Phys. Lett. 93, 212101 (2008)] M. Cson- 
tos, Y. Komijani, I. Shorabalko, K. Ensslin, D. Reuter, and A. D. Wieck, 
|Appl. Phys. Lett. 97, 022110 (2010) C. H. L. Quay, T. L. Hughes, J. A. 
Sulpizio, L. N. Pfeiffer, K. W. Baldwin, K. W.West, D. Goldhaber-Gordon, 
and R. de Picciotto, Nat. Phys. 6, 336 (2010). 

6 P. Y. Yu and M. Cardona, Fundamentals of Semiconductors, 2nd ed. 
(Springer, Berlin, 1999). 

7 J. M. Luttinger, Phys. Rev. 102, 1030 (1956). 

8 A. Baldereschi and N. O. Lipari. Phys. Rev. B 8, 2697 (1973). 

9 I. Vurgaftman, J. R. Meyer, and L. R. Ram-Mohan, J. Appl. Phys. 89, 5815 
(2001). 

10 We neglect the typically very small anisotropic contribution to Zeeman 

splitting in the valence band. 
"See, e.g., U. Zulicke, phys. stat. sol. (c) 3, 4354 (2006), where a similar 

model for hole point contacts was discussed. 
12 S. Datta, Electron Transport in Mesoscopic Systems (Cambridge U Press, 

Cambridge, UK, 1995). 
13 M. Biittiker, H. Thomas, and A. Pretre, Z. Phys. B 94, 133 (1994). 
14 P. W. Brouwer, [Phys. Rev. B 58, R10135 ( 1998) 

l5 These partial currents satisfy Onsager-type relations i t ( l )(B x ) = 
i t (,)(-B x ), it(By) = i t ,(-B y ), i r (,)(By) = i r(l) {-B y ), i r {B x ) = 



4 



V(— B x ), where the last two equations are valid for symmetric devices 
with the barrier placed midway between contacts. 
'N. R. Ogg, Proc. Phys. Soc. 89, 431 (1966). 



7 F. Malcher, G. Lommer, and U. Rossler, Superlat. Microstract. 2, 267 
(1986). 

8 We used the conduction-band-structure parameters given in M. Stadele, 
K. Hess, T. Ruf, and M. Cardona, J. Appl. Phys. 91, 9435 (2002). 



